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We present a new calculation of the CP violation parameter e'/e. 
The results reported in this paper have been obtained by using the 
AS = 1 effective Hamiltonian computed at the next-to-leading order, 
including QCD and QED penguins. The matrix elements of the rel- 
evant operators have been taken from lattice QCD, at a scale [i = 2 
GeV. At this relatively large scale, the perturbative matching between 
the relevant operators and the corresponding coefficients is quite reli- 
able. 

The effect of the next-to-leading corrections is to lower the predic- 
tion obtained at the leading order, thus favouring the experimental 
result of E731. We analyze different contributions to the final result 
and compare the leading and next-to-leading cases. 
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In this paper we present a theoretical prediction of e'/e obtained from 
the effective weak hamiltonian Ti^ft^, up to next-to-leading QCD and QED 
corrections. The Wilson coefficients of the operators of 'H^ff 1 have been 
computed using the (10 x 10) anomalous dimension matrix which governs the 
mixing of the relevant current-current and penguin operators, renormalized 
in the MS t'Hooft-Veltman dimensional regularization scheme (HV). The 
anomalous dimension matrix includes orders (a s t) n , a s (a s t) n , a e t(a s t) n and 
a e (a s t) n (where t = lnM^//i 2 ). The coefficients have been obtained by 
integrating numerically the renormalization group equations from a scale 
/i~M^y~m 4 down to a scale fi = 2 GeV > m c . The coefficients at the initial 
scale fi ~ My/ are those computed by Inami and Lim in ref . and by Flynn 
and Randall 0, as corrected in ref. ||. The initial effective hamiltonian 
includes QCD, Z° and electro-magnetic penguins and box diagrams. We 
have neglected the running of the coefficients between m t and M w , but this 
corresponds to a negligible error in the final result. 

The matrix elements of the operators, expressed in terms of the so-called 
5-factors, have been taken, when possible, from lattice calculations. It turns 
out that the most important 5-factors have indeed been determined on the 
lattice, with either Wilson and staggered fermions One notable ex- 

ception is the S-factor of the operator 0~ , believed to be responsable for the 
AI = 1/2 enhancement. For this reason we take the non-leptonic AI = 1/2 
amplitude from experiments. For those matrix elements not yet determined 
by lattice calculations, we allowed a variation of the i?-factor in the range 
1-6. The next-to-leading relation between the operators in the lattice renor- 
malization scheme and MS naive dimensional regularization (NDR) is known 
[fLCfl-[1l5| (the corrections are of the order of 30% in current lattice calcula- 
tions) and is usually included in the lattice results. We have only computed 
at one loop the factors which relate the operators renormalized in NDR to 
the corresponding ones in HV. These corrections give a negligible shift to the 
values of the B-parameters. Given the uncertainties in the evaluation of the 
B-parameters, we have decided to use the same values as in ref. |16[], see 
Table |l]. We have chosen the renormalization scale /i = 2 GeV, corresponding 
to the typical inverse lattice spacing at which actual calculations of weak 
matrix elements are performed in lattice QCD (a _1 = 2— 3.5 GeV). 

There are great advantages in using weak matrix elements of operators 
from lattice QCD. From the theoretical point of view, the matching of the 
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coefficients with the matrix elements of the operators, renormalized at the 
scale /i, is exact at the next-to-leading order in a s . This is not the case with 
other methods. For example, in the 1/N expansion, one has to match the 
coefficients, computed by renormalizing the operators on quark states, with 
the matrix elements calculated in the meson theory. However, it is not clear 
to us how this matching can be implemented at the next-to-leading order. 

A second important point is that the scale fi can be taken as large as 
2 — 3 GeV, where the perturbative evaluation of the coefficients is expected 
to be accurate and the final result quite stable for variations of fi in the above 
range, see below. This is to be contrasted with other approaches where values 
of \x ~ 0.6 — 0.8 GeV are chosen. At such low scales the results are not stable 
against a variation of \i. Moreover, at low values of /j,, next-to-leading order 
(NLO) corrections to the coefficients of penguin operators are large, so that 
one can question the convergence of the perturbative expansion. 

The calculation of e'/e has been combined with a next-to- leading order 
calculation of e and of the B°-B° mixing amplitude, following the approach of 



ref . []T6| . For these quantities the next-to- leading perturbative corrections to 
the Wilson coefficients have been known for quite a while and we do not have 
much to add to previous analyses. We notice that, with the inclusion of the 
coefficients and of the anomalous dimensions computed at second order, the 
analysis of ref. |I6J is here consistently done at the next-to-leading accuracy. 



The strategy is the same as in ref. |16] : from the comparison of the theoretical 



value of e with the experimental number, given the uncertainties on the 
matrix element of the AS = 2 operator and on the CKM parameters, we 
find a range of allowed values of the CP violating phase 5. Correspondingly 
we compute e'/ e, which will also be affected by a theoretical error, see Figs. 
(©-(§ below. 

The paper is organized as follows. We first give the effective Hamiltonians 
responsible for AS* = 2 and AS = 1 transitions, expressed in terms of the 
relevant operators and their Wilson coefficients. From these Hamiltonians we 
derive the expressions for e and e'/e, written as combinations of the Wilson 
coefficients times the B-parameters of the different operators, i.e. their matrix 
elements. We do not discuss the determination of the coefficients of the 
operators relevant for e, since this point was explained in great detail in ref. 
[frHl. We instead focus on the evaluation of the coefficients of the AS = 1 
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Hamiltonian beyond the leading order including the effects of the electro- 
weak penguins. A comparison of the predictions in the leading and next- 
to-leading cases is presented. We also study the uncertainties in the final 
evaluation of e'/e coming from the choice of the scale and Aq C d and the 
dependence on the top quark mass. Finally the theoretical predictions will 
be confronted with the experimental results coming from NA31 [|17]] and E731 



We have also considered the B -B transition parameter Xd, following the 
analysis done in ref. [16|, but with a different range of Aqq D = (340 ± 120) 



MeV |T9| , |2"0"| |. It remains true that a large value of for m t > 140-150 
GeV, favours a positive value for cos 5. Since the comparison of x^ with the 
theoretical prediction has no other effect on e'/e we will not discuss it any 
more. 

1) e : The effective Hamiltonian governing the AS* = 2 amplitude is given 

by: 

</ M = ^hK(dlZs) 2 {X 2 c F(x c ) + X 2 F(x t ) + 2X c X t F(x c ,x t )} 



(1) 



where Gp is the Fermi coupling constant and 7^ = 7^(1 — 75) ; A g 's are 
related to the CKM matrix elements by \ q = V^Vqf, where 'i' and '/' are 
the labels of the initial and final states respectively (in the present case i = s 
and / = d). x q = m 2 q / and the functions F(xi) and F(xi,Xj) are the 
so-called Inami-Lim functions 0, obtained from the calculation of the basic 
box-diagram and including QCD corrections [TL\. F(xt) is known at the 



next-to-leading order, which has been included in our calculation. From eq. 
([!]) we can derive the CP violation parameter e: 



\e\ e=0 = C e B K A 2 X G a {f(x c , x t ) + F{x t )[A 2 \\l - p)\ - F{x c )} (2) 



where 



= G\p K m K M 2 w 
6 QV2tt 2 AM k 1 ' 

AMk is the mass difference between the two neutral kaon mass eigenstates. 
In eq. (0) a — V P 2 + V 2 s ^ n $ an< i A, p, r] and S are the parameters of the CKM 
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matrix in the Wolfenstein parametrization P^| . Bk is the renormalization 
group invariant B-factor, defined as : 

B K = B K (p) [a QCD ^T 6/25 (4) 

Bk takes into account all the possible deviations from the vacuum inser- 
tion approximation in the evaluation of the (K°\(d'ji j s) 2 \K ) matrix element 
{Bk = 1 corresponding to an exact vacuum insertion approximation). 

2) e'/e: Most of the discussion in this paper is devoted to the Wilson 
coefficients of the operators appearing in the effective AS = 1 Hamiltonian, 
which we have computed at the next-to-leading order, including QCD and 
QED corrections. 

The AS* = 1 effective hamiltonian is given by: 

w$ =1| = £G(/*XM/*) (5) 

i 

The complete basis of operators when QCD and QED corrections are taken 
into account is given by: 



Qi = 


(s a d a t 


] (V-A)(Uf3Uf3)(V-A) 


Q 2 = 




(V-A)(upU a )(V-A) 


^3,5 = 


(s a d a t 


] (V-A) Y (Wp){vta) 






q=u,d,s 


Qlfi = 


(s a dp) 


(V-A) Y (.Wo){VtA) 
q=u,d,s 


Q7,9 = 


^(s a d 


a) (V-A) Y e q(.Wp)(V±A) 






q=u,d,s 


^8,10 — 


\is a d 


&){V-A) Y e q(Wa)(V±A) 






q=u,d,s 


Qi = 


(s a d a i 


] {V-A){Cf3Cp)(y-A) 


Q c 2 = 


{Sadp) 


(V-A)(c/3C a )(y-A) 



where the subscript (V ± A) indicates the chiral structure and a and f3 are 
colour indices. 

The operators Qi(fi) are renormalized at the scale \x in MS , using the 
HV regularization scheme. The corresponding coefficients, C;(/i) are scheme 
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dependent. The dependence on the regularization scheme appears at one 
loop, when we express the original current-current product in terms of the 
Wilson OPE: 

<F|^ =1| |/> = g 2 w /%jd i xD w {x\Ml l )<F\T(j^x) 1 Jl{0))\I> 
- J2 C i( M w)<F\Qi(M w )\I> (7) 

i 

The Wilson coefficients C{Mw) = (Ci(Mw), C^(Mw)...) are found by match- 
ing, at 0(a e ) and 0(a s ) in HV, the current- current and penguin diagrams 
computed with the W and top propagators to those computed with the local 
four-fermion operators in the effective theory. 

C(fi) are expressed in terms of C(Mw) through the renormahzation evo- 
lution matrix W[fi, Mw]fj[- 

d(ji) = W\n,M w ]C(M w ) (8) 

where: 

W[fi, M w \ = M[n]U[n, M W ]M'[M W ] (9) 

with: 

*M = ( i + S*)( 1 + ^)( 1 + ^) < 10 » 

and 

M'[M W ] = ( 1 - -^—P] ( 1 - (i - ?1k) (11) 



a s (M w ) J \ Arc J \ Att 

At the next-to- leading accuracy W[pt, Mw] is regularization scheme de- 
pendent. 

Eqs. (f|), (f[Qf ) and ([0]) require a detailed explanation. At the leading 
order, the QCD anomalous dimension matrix, including QCD penguins, has 



been computed in refs. [|23i |24j] . The electro- weak anomalous dimension 



matrix at the same order can be found in refs. H, [^5] and [p6 |. We have 



computed the anomalous dimension matrix at the next-to-leading order, by 



1 We have properly taken into account the beauty threshold in the evolution matrix. 
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calculating all the current-current and penguin operators at two loops up to 
order ot 2 s t and a e a s t. This corresponds to all the diagrams with four external 
quark legs, where one of the operators in the list given in eqs. (|5|) is inserted 
and two gluons or one gluon and one photon are exchanged. At 0(a^t), the 
two loop anomalous dimension matrix was computed in refs. |27], for 
current-current diagrams and in ref. |2{| for penguin diagrams. The explicit 
expression of the anomalous dimension matrix alone would take more space 
than that allowed for a letter and will be presented in a separate publication 
[poll . Here we simply explain the meaning of all the terms appearing in eq. 
(PD-(|TT|). To obtain the expression in eq. (|9|)-(|TT|) we have neglected the 
running of the coefficients between the top quark mass and the W mass. We 
have also expanded the formula at first order in a e and neglected the running 
of the electro-magnetic coupling. These approximations are immaterial for 
the final numerical result. 

The matrix U[fi, Mw] in eq. (|9|) is given by 



E>[//,M, 



a s {M ] 



w i 



(12) 



The matrices P, J and K are solutions of the equations: 

■MT- 



P + 



P 



Wo 



j(0)T 

~Wo~ 



(13) 



J 



J 



' 2/3 



#La(0)T _ 7. 



(1)T 



2 T« 



2f3o 



(14) 



01 



KM 0)T =li l)T + li 0)T J + li 1)T P+ li 0)T ,JP -2AP-£p7f r (15) 

The anomalous dimension matrix, which includes gluon and photon correc- 
tions has been separated in several pieces which appear in the above equa- 
tions: 
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^_s_ - (0) , <^e_ - (0) _|_ ( a s\2^(l) _|_ ^s_( 

47T 7s 47T 7e l 47T J 7s 47T47T 7 ' 



(1) 



(16) 



where each of the ■yf'^ is a 10 x 10 matrix. In eqs. ( |T2| - p~5D , (3$ and (5\ are 
the first two coefficients of the /3-function of a s . 
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From Ti)^ff 1 ' we can derive the expression for e'\ 

uj-\\mA 2 )' - (1 - fi JB ) ImAol (17) 



e 



ReA 

where (Imv4 2 )' and ImA are given by: 

lmA = -^/m(^l/ u ) {- (c 6 B 6 + \c 5 B 5 ) Z + (c 4 B 4 + \c 3 B 3 ) X+ 



y/2 "V" m J I V ° 3 
/2Y Z 



1/2 /2F Z X\ _ r1/2 /_ . Z . X\ 



^ 9 1/2 f + (18) 



C 8 5 8 3 / 2 (F-|) + C 9 5f^} (19) 

c<j = ReA 2 /ReA = 0.045 and we have introduced (Iim4 2 )' defined as: 

ImA 2 = (Irm4 2 )' + n IB (u ImAo) (20) 

Q IB = +0.25 ± 0.10 represents the isospin breaking contribution, see for 
example ref. |3Tfl . 

The numerical evaluation of e'/e requires the knowledge of the Wilson 
coefficients of the operators and of the corresponding matrix elements. The 
Wilson coefficients have been evaluated, using eq. (§), combined with the 
evolution matrix of eq. @, and the initial conditions computed in refs. |fj. § 
(and given for HV in ref. p9fl). The matrix elements of the operators have 
been written in terms of the three quantities (see eqs. fll8D and ( |T9|) ) : 

X = U(M 2 K -Ml) (21) 



and 



M 2 K \ __/0.15GeV' 



2 



Y = f* mi n ~UX[ - (22) 
\m 8 {fj,) + m d {fj,) J \ m s (fx) J 

Z = 4 (y- ( 23 ) 
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r d(3/2) rc d d r (1/2) r (3/2) 



0.8 ±0.2 0-0.15W 1-6W 1.0 ±0.2 1« 1.0 ±0.2 



Table 1: Values of the 5-parameters. Entries with a M are educated guesses; 
the others are taken from lattice QCD calculations. 



and a set {-Bj} of B-parameters (in our normalization f n = 132 MeV). The 
numerical value of the B-parameters have been taken from lattice calculations 
and multiplied by suitable renormalization factors to take into account the 
difference between HV and the lattice regularization scheme. For those B- 
factors which have not been computed yet on the lattice we have used an 
educated guess, which will be discussed below. We observe that in eqs. ([1S|) 
and ( |i9|) only nine coefficients (B-parameters) appear since we have used the 
relation Q 10 = —Q 3 + Q 4 + Q 9 . 

We will call "central" results obtained by using the central values of 
the B-parameters reported in Table |], fi — 2 GeV and Aqcd = 340 MeV 
(4 flavours). We will study the uncertainty of the theoretical prediction by 
varying the B-parameters, m s and the experimental quantities, such as Aqcd, 
the CKM mixing parameters A, p, etc., in the range indicated by the errors 
on the quantities reported in Tables 1 and 2. In the study of the differences 
between e'/e, computed at the leading order (LO) or next-to-leading order 
(NLO), we will keep fixed the coefficients at the scale My/. This means that, 
also at the LO, the coefficients include the 0(a s ) and 0(a e ) corrections which 
arise from the matching between the effective hamiltonian and the original 
current-current product. We also keep fixed Aq C d m passing from the LO 
calculation to the NLO one. 

We first consider the relative contribution of different operators to e'/e. 
Following the standard notation we write: 

e'/e~RxC 6 B 6 (l-J2^) (24) 
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Figure 1: Band of allowed values for e'/e at m t = 100 GeV (next-to- 
leading order). The dashed lines represent the experimental results of NA31, 
(2.3 ± 0.7) ■ 10~ 3 and E731, (0.6 ± 0.7) • lO" 3 . 



where the coefficient R includes sin 5. In eq.(24) fij = fi^fi 3 / 2 , etc. indicate 
the relative corrections due to the matrix elements of the operators, different 
from Q 6 , appearing in eqs. (|18D-(|19|). For example, the central value of R, 
averaged over the values of cos S allowed by the analysis of e, is \R\ ~ 1.0- 10~ 2 
for m t = 140 GeV. 

In the following discussion we fix /i = 2 GeV and m t = 140 GeV and 
we allow a variation of the B-parameters, computed on the lattice, around 
their central values. For those matrix elements still to be computed we 

1 /2 

proceed as follows. We fix B 7 _ 9 equal to one, since their contribution is of 
the order of few percent. If we take B^^ = 1 and B\ 2 = 0, at the LO the 

largest contributions to e'/e come from Q$, Q 8 and Q$, £l^ 2 ~ 28% and 
Slg^ 2 ~ —18%. Cg is much larger than Cg and compensates for the matrix 
element of Q$ which is much smaller than the matrix element of Q$. With 
i?3 5 4 = 1, ^3 ~ —0.01 and ~ 0.04. The operators Q^^ have a chiral 
structure similar to Q2 and, at scales larger than m c , give rise to the same 
"eye" diagrams, like Q2 does. In order to explain the experimental AI =1/2 
enhancement, B 2 must be of the order 5-6. For this reason in ref. [|16j and 
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Figure 2: Same as in Fig. (jT|) at m t = 140 GeV. 



in the present calculation we have allowed a variation of -83,4 between 1 and 
6. The central values in this case are Q% ~ —0.04 and Q4 ~ 0.28. Thus, if 
Q4 has a large B-factor, its contribution is as important as the contribution 
of Qs or Q 9 . B\ 2 are zero in the vacuum saturation approximation. Their 
coefficients however are so large that even a small B-factor can give a sizeable 
contribution to e'/e. By varying Bf 2 between 0-0.15 we find that fl 2 ~ —0.18 
(fll ~ 0.02), comparable to other large terms. 

In the following we will focus on the main contributions to e'/e, which 
are due to few operators, (54,6,7,8,9 and Q 2 . All the other terms are of the 
order of a few %, with alternating signs and we will not discuss them any 
more. 

We now come to the effects of the next-to-leading corrections. If we only 
include next-to-leading corrections due to two gluon exchanges (correspond- 
ing to 7^), all the coefficients are slightly changed, without modifying the 
pattern of the different contributions observed at the leading order. For ex- 
ample fi 4 goes from 0.28 to 0.34, fi 5 from 0.06 to 0.08 and fig /2 from 0.28 to 
0.26. Vt^ 2 increases from —0.18 to —0.22. The gluon-photon corrections on 
Qi are always very tiny and, given the uncertainties on the B-factors, practi- 
cally invisible. The only exceptions 7 is very small at the leading 
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Figure 3: Same as in Fig. (jT|) at m t = 200 GeV. 



order, —0.01, becomes larger including two-gluon exchanges, —0.05 and be- 
comes —0.09 with the complete next-to-leading corrections; f2 8 is 0.28 at 
the LO and increases to 0.51 with the complete NLO corrections. We notice 
that the sum of the contributions of Q^^ and Qg almost cancel at the lead- 
ing, as well as at the next to leading order. f2^ 2 + Q^ 2 + Q^ 2 ~ 0.10 at the 
leading order and 0.20 at the next-to-leading one. Similarly Q4 + Q2 ~ +10% 
both at the LO and at NLO. We thus find that at m t = 140 GeV and for the 
"central" values of the B-parameters, the result is essentially the same as the 
original Gilman-Wise prediction, since the most important electro-penguin 
corrections almost cancel and their sum is of the order of a few per centQ. 
Thus e'/ e is essentially determined by Cq and the corresponding B-parameter, 
B 6 . C 6 varies from -6.2 x 10~ 2 to -5.1 x 10~ 2 , at fixed A QCD = 340 MeV, so 
that the average value of e'/e decreases from 6.5 x 10 -4 ( LO) to 4.5 x 10~ 4 
(NLO). The decrease of Cq is mainly due to the variation of the value of a s 
between LO and NLO, if we insist in using the same value of Aqcd- The rel- 
ative uncertainty on Cq, by using different values of Aqcd, remains roughly 

2 The cancellation between different contributions may not occur for values of the B- 
parameters different from their central ones. This will correspond to a band of uncertainty 
in the theoretical prediction for e'/e, which is reported in the Figs. (Q)-(|4|). 
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Figure 4: Same as in Fig. (|T[) at m t = 100 GeV, with the coefficients of the 
operators for e'/e computed at the leading order. 



the same between the LO, 5C 6 /C 6 ~ 18% and NLO, 5C 6 /C 6 ~ 17%. 

As observed already in ref. 0, and confirmed by all other analyses 



16], for increasing values of the top mass, the contribution of the electro- 



penguin operators tends to cancel the contribution coming from Q e . We 
have already observed that at m t = 140 GeV the main corrections come 
from fi^g 9, ^4 and f2|. f2 4 and Vt^ vary very little with m t and their values 
remain essentially the same in going from LO to NLO. We know very little 
about their matrix elements from non-perturbative calculations. We estimate 
that the uncertainty on the final result, coming from the poor knowledge of 
the matrix elements of and Q\ is of the order of 30%. On the other hand 
^7 8 9 vary with m t and change considerably from LO to NLO. For example, 

at m t = 140 GeV Slf 2 is ~ 0.28 at LO and ~ 0.50 at NLO, at m t = 200 GeV 
it becomes 0.68 at LO and 0.90 at NLO. We also observe large variations for 
n 3 7 /2 and fig /2 Q. The sum tt 3 7 /2 + tf 8 /2 + fig /2 is ~ 20% at m t = 140 GeV 
and ~ 65% at m t = 200 GeV. The same can be said also at low values of 
the top mass, for example m t = 100 GeV. In summary, in the range of mt 



3 The increase of the relative contribution of fi^g is due to the decrease of Cq combined 
with an increase of C7_g between the leading and next-to-leading cases. 
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parameter 



value 



A-qcd 
m s (2GeV) 
m c (2 GeV) 
m b {2 GeV) 
A\ 2 

Vp r +Tf = v ub /(\v cb ) 

^exp 

ReA 



340 ± 120 GeV 
(170 ± 30) MeV 
1.5 GeV 
4.5 GeV 
0.047 ±0.004 
0.50 ±0.14 
2.28 • 10~ 3 
2.7- 10~ 7 GeV 



Table 2: Values of experimental parameters used in this work. 



considered in this work, even though single contributions from the electro- 
penguin operators may change by ~ 40% because of the next-to-leading 
corrections, they give globally more or less the same relative contribution at 
NLO as at the leading order. Since for fixed Aqcd, Cq decreases, the net effect 
is the the central value of the theoretical prediction for e'/e is smaller at the 
NLO. In Figs. (|l])-(|3|), we report our results at NLO for e'/e, at three different 
values of m t , 100, 140 and 200 GeV. The predictions are given for different 
values of cos 5, which are compatible with the analysis of e [16]. They have 
been obtained by comparing the theoretical expressions for e, eqs. (@), to the 
corresponding experimental value. The dashed lines indicate the results by 
the NA31 and E731 experiments at the 1-a level. The numbers and errors 
reported in the figures are the average and the theoretical variance (for cos S 
positive or negative) computed by varying the B-parameters and the other 
quantities in the ranges reported in Table [TJ and |2|. Since the value of e'/ e is 
lowered at the next-to-leading order, our theoretical prediction, at m t = 140 
GeV, is now centered on the experimental result of E731 (cos 5 > 0)0- For 
a comparison between the LO and the NLO, we report in Fig. (|]) the LO 
result, at mt = 100 GeV. The LO result in this case sits in the middle between 



4 We recall that the most recent lattice results 32 -|34] and QCD sum rules calculations 




suggest a large value of Jb which corresponds to cosS > when m t > 140-150 

GeV ^ 
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Figure 5: C 6 and C$ as a function of /i for Aq C d = 220 (dotted), 340 (solid) 
and 460 (dashed) MeV. 

the measurements of NA31 and E7310. We also observe that the band of 
error is slightly reduced by the inclusion of the next-to-leading corrections. 
Indeed the relative error on e'/e at LO and NLO is basically the same. Since 
the central value of e'/e is decreased at NLO, the band in Fig. ([[]) appears 
narrower than in Fig. (f|). 

Before concluding this paper we want to discuss the stability of the Wil- 
son coefficients, i.e. how much their values are affected by the presence of 
NLO corrections and their sensitivity with respect to a variation of and 
Aqcd- We observe that for ji < 1 GeV the values of the coefficients start to 
vary wildly if one changes Aq C d or /i. As an example, in Fig. (|5]), we report 
the variation of C 6 (/z) and C 8 (/i) as a function of /x, at three different values 
of Aqcd and mt = 140 GeV. The coefficients change by a factor of 2 or more, 
for fi < 1 GeV, if we vary Aqcd from 220 to 460 MeV. Moreover, at small 
values of /i, the difference between the leading and next-to-leading results 
is very large, and we cannot trust the perturbative expansion. To ilustrate 
this point, we report in Fig. @ C 8 (fi), for 0.8 GeV < /i < 4.0 GeV at the 
leading and next-to-leading order. On the basis of the above discussion we 

5 The corresponding figure in ref. (l6) was obtained with different values of Aqcd- 
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Figure 6: C 8 as a function of /i for A QCD = 340 MeV at the LO and NLO. 

believe that a realistic prediction for e'/e can only be obtained by matching 
the Wilson coefficients and the operators matrix elements at scales larger 
than m c , as one can do by taking the matrix elements from lattice QCD. 

The conclusions of this work are the following. It is reassuring that by 
taking a renormalization scale \x of the order of 2 GeV, the Wilson coefficients 
of the effective Hamiltonian do not vary by more than ~ 35% and that the 
final result on e'/e is quite stable in going from the leading to the next-to- 
leading order. At NLO one gives a precise meaning to the value of Aq C d to 
be used. It is now consistent to take Aqqd from the measurements done in 
deep inelastic scattering or at LEP. The NLO calculation of the anomalous 
dimension matrix makes also consistent the evolution of the Wilson coeffi- 
cients with their matching at /i ~ m t , Mw- Indeed the matching procedure 
is a next-to-leading order effect. Finally, in the operator product expansion, 
only the NLO calculation fixes unambigously the scale at which the operators 
must be renormalized. In the lattice case the scale is dictated, up to higher 
order effects, by the inverse lattice spacing at which the numerical simula- 
tions are performed. The major source of theoretical uncertainty remains 
now the evaluation of the matrix elements, which are determined with large 
errors or are still to be determined, like it is the case for 0~. 
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From the phenomenological point of view, at fixed Aqcd and nit, the 
next-to-leading corrections lower the theoretical predictions, thus favouring 
the experimental result by E731. 
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